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The absence of compelling theoretical model requires the parameterizing the dark energy to probe 
its properties. The parametrization of the equation of state of the dark energy is a common method. 
We explore the theoretical optimization of the parametrization based on the Fisher information 
matrix. As a suitable parametrization, it should be stable at high redshift and should produce 
the determinant of the Fisher matrix as large as possible. For the illustration, we propose one 
parametrization which can satisfy both criteria. By using the proper parametrization, we can 
improve the constraints on the dark energy even for the same data. We also show the weakness of 
the so-called principal component analysis method. 

PACS numbers: 95.36.+X, 95.80.+p, 98.62.Ve. 

We consider strategies for the most accurate determination of the dark energy (DE) parameters using 
cosmological observables. Due to the absence of compelling theoretical model, the parameterizing the DE 
by its equation of state (eos) u> is commonly used in the analysis. It is usually considered what is the optimal 
redshift distribution to best constrain those parameters [1, 2]. However, we contemplate the theoretical 
optimization by using the Fisher information matrix and compare the different parametrization of uj. 

If one has a dataset of observable Oi with redshifts Zi, i = 1, ■ • • , N, one may compute the least squares 
estimates of the cosmological parameters by minimizing 



Xo = 



E- — -2 — (o- 1 ) 



where Oi is the measured quantity, 0(zi,p) is the expected value of the observable for the redshift Zi, p is 
the set of cosmological parameters to estimate, and Oj is the error on the measurement. Then, the Fisher 
information matrix is defined to be 
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where p* m are the parameter estimates (i.e. where Xo becomes minimum). Since the likelihood function 
is approximately Gaussian near the maximum likelihood (ML) point, the covariance matrix for a maximum 
likelihood estimator is given by 
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Pi* >Pm 



The Fisher information matrix is simply the expectation value of the inverse of the covariance matrix at 
the ML point. From now on, we will omit the square bracket and the evaluating parameter estimates in 
the Fisher matrix for the convenience. 

The iso-Axo contour in the parameter space is approximated by the quadratic equation 

(Ap) T FAj?= A x 2 o, (0.4) 
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where Ap = p — p* is the deviation of the parameters from the fiducial value and F is the Fisher matrix 
given in Eq. (0.2). Thus, the most accurate determination of n cosmological parameters is described by a 
minimal volume of the n-dimensional ellipsoid. The volume of the ellipsoid is given by 
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Thus, the determination of F should be maximum to minimize the volume of the uncertainty ellipsoid 
By using Eqs. (0.1) and (0.2), we obtain 
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where we use ^ — — ^j z "^ ^3^3^'^ — and define the approximate Fisher matrix Fi m including only the 
second term in the original i*j m . We can use this approximation because \ 2 becomes minimum at ML point 
P 
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Eq. (0.7) can be interpreted as the summation of 



weighted by 



— - — i ' is zero. As long 



as the data converges to the fiducial values, the weighted summation of the variation of 



dO(zj,p) 



(i.e. 



^2 °* °Jf l 'V ) d Qp i l Qp ) should be negligible. Thus, the first term of the right hand side of Eq. (0.6) is much 
smaller than the second one. We will explicitly show this later in Table. II by using H(z) data. 

In general, p includes the present value of the Hubble parameter Hq, the present value of the matter 
energy density contrast O m o, the present value of the DE density contrast ^deo and the parameters of the 
eos of DE to. We concentrate on the flat universe and firjEo i s determined as 1 — fi m0 . Now, we assume that 
the nuisance parameters Hq and f2 m o are exactly known and only consider the estimation of the parameters 
of u. F given in Eq. (0.6) can be further specified depending on the observables. First, we consider when 
Pde is a direct variable of 0( Z i,p). H(z) corresponds to this case. Then, the approximate Fisher matrix 
becomes 



Fir, 



dO(zj,p) 
^ af y d\n[p BE ( Zi ,p)} 



N 




— dln(l + x 

dpi 



dpm 



dln(l + x) 




dln(l + x) 



(0.8) 



where we use = pdeo ex P 
like in d^, then F is given by 
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Second, if pde is included in O as an integrand 
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0.09 


0.17 


0.27 


0.4 


0.48 


0.88 


0.9 


1.3 


1.43 


1.53 


1.75 


H(z) 


69 ± 12 


83 ±8 


77 ± 14 


95 ±17 


97 ±62 


90 ±40 


117 ± 23 


168 ± 17 


177 ± 18 


140 ± 14 


202 ± 40 



TABLE I: H(z) from passively evolving galaxies. 



In general, we can express the approximate Fisher matrix as 

N n n 



Gi = 



E E E GlWi{zi)W m {zi) , where 

i— 1 l — l m—1 

( Hp 3[(g( gi )/go) 2 -n m o(l+ Zi ) 3 ] 
V 



2(H(z,)/H ) 



if O = H{z) 



and 



Wi{z t ) 



Jo 



dln(l ± x) 



l-Zi 30O(x) 



\iO = d L {z) 

if O = H{z) 
f*MMl d \n(l + y)\dx if = d L (z). 



(0.10) 



r-Zi it) 
L JO d\n[p DE (x)] 

Thus, the determinant of the approximate Fisher matrix is given by 
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where e a 6c- -j analogous to the Levi-Civita symbol, is ±1 for even commutations (abc- ■ ■ ) and —1 for odd 
commutations, and zero if any index is repeated. 

First, we consider the models with two parameters uj = uj a ± ui b f{z). In this case, det(-F) is given by 
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If we choose the observable as H(z), then det(F) becomes 
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We use the -ff(z) data obtained from the passively evolving galaxies [3]. This is shown in the table I. 



In Eq. (0.14), we consider three well known parametrization [4-7]. This equation shows the importance of 
the parametrization of DE eos. The values of Gi and Gj should be very close to each other for the different 
models and thus the magnitude of det(i ;l ) depends only on the geometric part of w. Thus, the magnitude 
of det(-F) depends on the magnitude of the quantities in the square bracket. It means that we can extract 
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{a, a) 


(a,b) 


(6,6) 


(a, a) 


(a,b) 


(6,6) 


z 


0.19 


-0.66 


-1.04 


92.51 


52.19 


31.29 


-0.834 


0.211 


9.486 


149.3 


[4] 


ln[l + z] 


-0.12 


-0.46 


-0.45 


93.14 


38.52 


16.72 


-0.875 


0.388 


9.446 


65.07 


[ ] 


z 

1 + z 


-0.34 


-0.33 


-0.21 


93.66 


29.65 


9.75 


-0.916 


0.638 


9.412 


31.19 


[6, 7] 


tanh[ln[l + z]] 


-0.34 


-0.41 


-0.31 


93.63 


34.66 


13.38 


-0.906 


0.520 


9.413 


46.35 





TABLE II: Comparison of analysis results for the linear, logarithmic, CPL, and tanh[ln[l + z]] parametrizations by 
using H(z) data given in Ref. [3]. The components of the first term of Fisher matrix is negligible compared to the 
second one for all the models, to* and ujI are the best fit values of u) a and uib, respectively. Xmin are almost identical 
to each other. We use Hq = 71 km/Mpc/sec and f2 m o = 0.24. 



the better constraints depending on the paramctrization of u; even for the same measurement. f(z) = z or 
ln[l + z\ gives the stronger constraints on ui a and ui b compared to f{z) — However, both cases suffer 
the divergence of the error estimate as z increases. Given the covariance matrix, the error on lo is given by 
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C bb f(z) 2 + 2C ab f(z), 



where we use the two parameters models in the second equality. Thus, even though both the linear and 
the logarithmic parametrization of ui give the stronger constraints on parameters ui a and ui b as shown in 
the above, both approaches suffer from the instability of the error estimate as z increases. We propose the 
new parametrization of ui = ui a + Wfctanh[ln[l + z]] = uj a + uj b jjffij^p which is free from the divergence. 
Then the determinant of the approximate Fisher matrix is given by 
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This gives the larger det(-F) compared to CPL. 

We show the analysis results for the different parametrizations in table II. We use Hq = 71 km/Mpc/sec 
and fi m rj = 0.24 in the analysis. First, we confirm that ^ Oi-o^^v) d g ®(^J^ is negligible compared to 

E-4 d °( z "P> a °( z "P) m a n f models. We also show the best fit values of w„ and oj b for the different 

parametrizations. The values of the minimum \ 2 f° r the different models are almost same. But the values 
of determinant of Fisher matrix for the different models are quite different from one another. As we show 
in Eq. (0.14), linear parametrization produces the largest dct(F) and the so-called CPL parametrization 
(i.e. f(z) = yqjj) gives the smallest one among the models. 



Again, we show the above results in Fig. (1). In each column of the first row, we show the evolution of to 
for the different paramctrizatons. The thick solid lines are the best fit w and the shadow contour represents 
the 1-cr confidence level around the best fit. The linear and the logarithmic parametrizations diverge as z 
increases. Both CPL and tanh[ln[l + z}] are free from this divergence and both are almost identical. In 
the second row, we show the 1-cr (inner ellipsoid) and 2-cr (outer one) confidence contours in the oj a -u) b 
plane for the corresponding parametrizations. As we explain, the linear parametrization gives the strongest 
constraints on u a and tu b and the CPL produces the weakest constraints among the models. 

We repeat the analysis by using the so-called the principal component analysis (PCA). For the given H(z) 
data, we obtain ui a = —0.672, ui b = —0.783, ui c = —0.511, and u)d = —0.700 if we choose the 4 principal 
components at z a = 0.1, z b = 0.4, z c = 0.8, and Zd — 1.25. We emphasize that is induced values from 
the other three c^s as expalined in Ref. [8, 9]. If one does not use the constraints on ujj, then H(z) shows 
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FIG. 1: a) The variation of the uj over redshift for the different parametrizations. The thick solid lines are the 
best-fit and the shadow contour represents the 1-a confidence level around the best fit. b) The 1-a (inner ellipsoid) 
and 2-a (outer one) confidence contours in the ui a -u)b plane for the corresponding parametrizations to a). 
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FIG. 2: a) The evolution behavior of to for the different parametrizations b) H(z) obtained from the different 
parametrizations c) Relative errors of each H(z) from one obtained from u> a + uit tanh[ln[l + z]]. 

the discrepancy from the true value. We compare the evolutions of uj and H{z) obtained from each model 
in Fig. 2. In the first column, we show the evolution behavior of uj for linear (solid), logarithmic (dotted), 
CPL (dashed), tanh (dot-dashed), and PCA (long dashed), respectively. Both CPL and tanh shows the 
almost identical behavior for uj. We also show the evolution of H(z) including the data points in the second 
column. All of H(z) are too close to be distinguished. Thus, we compare the H{z) for the different models 
with one from ACDM model in the third column. The best fit values of H and il m0 in ACDM are 70.8 
km/Mpc/sec and 0.312, respectively with Xmm = 9.52. We use Hq = 71 km/Mpc/sec and fl m o = 0.24 in 
other two parameters models. AH means the difference of H(z) obtained from the parametrization and 
the ACDM. Except PCA, errors are smaller than 2 %. The fit uses the data up to z = 1.75 and that is 
why there exists bigger error at z > 2 for each model. 
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We can extend the previous consideration to models with n parameters and N datapoint. Then, the 
determinant of the approximate Fisher matrix is obtained from Eq. (0.11) 

N-n+lN-n+2 N / \ 2 

det(F)= S ■••X) G ? G ?"" G ?( s i3-iW a {zi)W b (z j )---W n (z l )\ . (0.17) 

i— 1 l—n \ / 

Thus, the determinant of the F is composed of (^) number of components with 

• • ■ G\ ^£ij—iW a (zi)Wb(zj) ■ ■ ■ W n {zi)\ where the binomial coefficient is given by (^) = 

N(N — 1)(N — 2) • ■ • (N — n + l)/n\. Each component has the nl number of permutations. If we 
consider the observable H(z), then the number of components inside the bracket of PCA method is smaller 
than n! depending on the location of the principal components and the rcdshift distribution of the data. 
It is due to the fact that W a {z) where a = a, is zero if data point is out of the corresponding principal 
component bin. Thus, PCA method gives less constraints to estimates than other paramctrization for 
the same number of parameters. We find two important features from this. First, PCA method has the 
weaker constraints than other methods because it produces smaller dot(F)u- Second, det(F) is sensitive 
to both the location of the PCA components and the data distribution . It is partly agree to the argument 
in Ref. [10]. 

We show the importance of paramctrization to have the better constrain on the dark energy. Even for 
the same data, we can improve the constraints of the dark energy if we find the proper parametrization of 
it. One can speculate the better parametrization based on this approach. Even though PCA is one of the 
most model independent parametrization, it suffer from the weaker constrain power than other methods. 
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